\(\int \frac {\cos ^3(c+d x)}{(e+f x)^2 (a+a \sin (c+d x))} \, dx\) [268]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [C] (verified)
   Fricas [A] (verification not implemented)
   Sympy [F(-1)]
   Maxima [C] (verification not implemented)
   Giac [C] (verification not implemented)
   Mupad [F(-1)]

Optimal result

Integrand size = 28, antiderivative size = 175 \[ \int \frac {\cos ^3(c+d x)}{(e+f x)^2 (a+a \sin (c+d x))} \, dx=-\frac {\cos (c+d x)}{a f (e+f x)}-\frac {d \cos \left (2 c-\frac {2 d e}{f}\right ) \operatorname {CosIntegral}\left (\frac {2 d e}{f}+2 d x\right )}{a f^2}-\frac {d \operatorname {CosIntegral}\left (\frac {d e}{f}+d x\right ) \sin \left (c-\frac {d e}{f}\right )}{a f^2}+\frac {\sin (2 c+2 d x)}{2 a f (e+f x)}-\frac {d \cos \left (c-\frac {d e}{f}\right ) \text {Si}\left (\frac {d e}{f}+d x\right )}{a f^2}+\frac {d \sin \left (2 c-\frac {2 d e}{f}\right ) \text {Si}\left (\frac {2 d e}{f}+2 d x\right )}{a f^2} \]

[Out]

-d*Ci(2*d*e/f+2*d*x)*cos(2*c-2*d*e/f)/a/f^2-cos(d*x+c)/a/f/(f*x+e)-d*cos(c-d*e/f)*Si(d*e/f+d*x)/a/f^2+d*Si(2*d
*e/f+2*d*x)*sin(2*c-2*d*e/f)/a/f^2-d*Ci(d*e/f+d*x)*sin(c-d*e/f)/a/f^2+1/2*sin(2*d*x+2*c)/a/f/(f*x+e)

Rubi [A] (verified)

Time = 0.23 (sec) , antiderivative size = 175, normalized size of antiderivative = 1.00, number of steps used = 11, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.250, Rules used = {4619, 3378, 3384, 3380, 3383, 4491, 12} \[ \int \frac {\cos ^3(c+d x)}{(e+f x)^2 (a+a \sin (c+d x))} \, dx=-\frac {d \sin \left (c-\frac {d e}{f}\right ) \operatorname {CosIntegral}\left (\frac {d e}{f}+d x\right )}{a f^2}-\frac {d \cos \left (2 c-\frac {2 d e}{f}\right ) \operatorname {CosIntegral}\left (\frac {2 d e}{f}+2 d x\right )}{a f^2}+\frac {d \sin \left (2 c-\frac {2 d e}{f}\right ) \text {Si}\left (\frac {2 d e}{f}+2 d x\right )}{a f^2}-\frac {d \cos \left (c-\frac {d e}{f}\right ) \text {Si}\left (\frac {d e}{f}+d x\right )}{a f^2}+\frac {\sin (2 c+2 d x)}{2 a f (e+f x)}-\frac {\cos (c+d x)}{a f (e+f x)} \]

[In]

Int[Cos[c + d*x]^3/((e + f*x)^2*(a + a*Sin[c + d*x])),x]

[Out]

-(Cos[c + d*x]/(a*f*(e + f*x))) - (d*Cos[2*c - (2*d*e)/f]*CosIntegral[(2*d*e)/f + 2*d*x])/(a*f^2) - (d*CosInte
gral[(d*e)/f + d*x]*Sin[c - (d*e)/f])/(a*f^2) + Sin[2*c + 2*d*x]/(2*a*f*(e + f*x)) - (d*Cos[c - (d*e)/f]*SinIn
tegral[(d*e)/f + d*x])/(a*f^2) + (d*Sin[2*c - (2*d*e)/f]*SinIntegral[(2*d*e)/f + 2*d*x])/(a*f^2)

Rule 12

Int[(a_)*(u_), x_Symbol] :> Dist[a, Int[u, x], x] /; FreeQ[a, x] &&  !MatchQ[u, (b_)*(v_) /; FreeQ[b, x]]

Rule 3378

Int[((c_.) + (d_.)*(x_))^(m_)*sin[(e_.) + (f_.)*(x_)], x_Symbol] :> Simp[(c + d*x)^(m + 1)*(Sin[e + f*x]/(d*(m
 + 1))), x] - Dist[f/(d*(m + 1)), Int[(c + d*x)^(m + 1)*Cos[e + f*x], x], x] /; FreeQ[{c, d, e, f}, x] && LtQ[
m, -1]

Rule 3380

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[SinIntegral[e + f*x]/d, x] /; FreeQ[{c, d,
 e, f}, x] && EqQ[d*e - c*f, 0]

Rule 3383

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[CosIntegral[e - Pi/2 + f*x]/d, x] /; FreeQ
[{c, d, e, f}, x] && EqQ[d*(e - Pi/2) - c*f, 0]

Rule 3384

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Dist[Cos[(d*e - c*f)/d], Int[Sin[c*(f/d) + f*x]
/(c + d*x), x], x] + Dist[Sin[(d*e - c*f)/d], Int[Cos[c*(f/d) + f*x]/(c + d*x), x], x] /; FreeQ[{c, d, e, f},
x] && NeQ[d*e - c*f, 0]

Rule 4491

Int[Cos[(a_.) + (b_.)*(x_)]^(p_.)*((c_.) + (d_.)*(x_))^(m_.)*Sin[(a_.) + (b_.)*(x_)]^(n_.), x_Symbol] :> Int[E
xpandTrigReduce[(c + d*x)^m, Sin[a + b*x]^n*Cos[a + b*x]^p, x], x] /; FreeQ[{a, b, c, d, m}, x] && IGtQ[n, 0]
&& IGtQ[p, 0]

Rule 4619

Int[(Cos[(c_.) + (d_.)*(x_)]^(n_)*((e_.) + (f_.)*(x_))^(m_.))/((a_) + (b_.)*Sin[(c_.) + (d_.)*(x_)]), x_Symbol
] :> Dist[1/a, Int[(e + f*x)^m*Cos[c + d*x]^(n - 2), x], x] - Dist[1/b, Int[(e + f*x)^m*Cos[c + d*x]^(n - 2)*S
in[c + d*x], x], x] /; FreeQ[{a, b, c, d, e, f, m}, x] && IGtQ[n, 1] && EqQ[a^2 - b^2, 0]

Rubi steps \begin{align*} \text {integral}& = \frac {\int \frac {\cos (c+d x)}{(e+f x)^2} \, dx}{a}-\frac {\int \frac {\cos (c+d x) \sin (c+d x)}{(e+f x)^2} \, dx}{a} \\ & = -\frac {\cos (c+d x)}{a f (e+f x)}-\frac {\int \frac {\sin (2 c+2 d x)}{2 (e+f x)^2} \, dx}{a}-\frac {d \int \frac {\sin (c+d x)}{e+f x} \, dx}{a f} \\ & = -\frac {\cos (c+d x)}{a f (e+f x)}-\frac {\int \frac {\sin (2 c+2 d x)}{(e+f x)^2} \, dx}{2 a}-\frac {\left (d \cos \left (c-\frac {d e}{f}\right )\right ) \int \frac {\sin \left (\frac {d e}{f}+d x\right )}{e+f x} \, dx}{a f}-\frac {\left (d \sin \left (c-\frac {d e}{f}\right )\right ) \int \frac {\cos \left (\frac {d e}{f}+d x\right )}{e+f x} \, dx}{a f} \\ & = -\frac {\cos (c+d x)}{a f (e+f x)}-\frac {d \operatorname {CosIntegral}\left (\frac {d e}{f}+d x\right ) \sin \left (c-\frac {d e}{f}\right )}{a f^2}+\frac {\sin (2 c+2 d x)}{2 a f (e+f x)}-\frac {d \cos \left (c-\frac {d e}{f}\right ) \text {Si}\left (\frac {d e}{f}+d x\right )}{a f^2}-\frac {d \int \frac {\cos (2 c+2 d x)}{e+f x} \, dx}{a f} \\ & = -\frac {\cos (c+d x)}{a f (e+f x)}-\frac {d \operatorname {CosIntegral}\left (\frac {d e}{f}+d x\right ) \sin \left (c-\frac {d e}{f}\right )}{a f^2}+\frac {\sin (2 c+2 d x)}{2 a f (e+f x)}-\frac {d \cos \left (c-\frac {d e}{f}\right ) \text {Si}\left (\frac {d e}{f}+d x\right )}{a f^2}-\frac {\left (d \cos \left (2 c-\frac {2 d e}{f}\right )\right ) \int \frac {\cos \left (\frac {2 d e}{f}+2 d x\right )}{e+f x} \, dx}{a f}+\frac {\left (d \sin \left (2 c-\frac {2 d e}{f}\right )\right ) \int \frac {\sin \left (\frac {2 d e}{f}+2 d x\right )}{e+f x} \, dx}{a f} \\ & = -\frac {\cos (c+d x)}{a f (e+f x)}-\frac {d \cos \left (2 c-\frac {2 d e}{f}\right ) \operatorname {CosIntegral}\left (\frac {2 d e}{f}+2 d x\right )}{a f^2}-\frac {d \operatorname {CosIntegral}\left (\frac {d e}{f}+d x\right ) \sin \left (c-\frac {d e}{f}\right )}{a f^2}+\frac {\sin (2 c+2 d x)}{2 a f (e+f x)}-\frac {d \cos \left (c-\frac {d e}{f}\right ) \text {Si}\left (\frac {d e}{f}+d x\right )}{a f^2}+\frac {d \sin \left (2 c-\frac {2 d e}{f}\right ) \text {Si}\left (\frac {2 d e}{f}+2 d x\right )}{a f^2} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.50 (sec) , antiderivative size = 203, normalized size of antiderivative = 1.16 \[ \int \frac {\cos ^3(c+d x)}{(e+f x)^2 (a+a \sin (c+d x))} \, dx=\frac {-2 f \cos (c+d x)-2 d (e+f x) \cos \left (2 c-\frac {2 d e}{f}\right ) \operatorname {CosIntegral}\left (\frac {2 d (e+f x)}{f}\right )-2 d (e+f x) \operatorname {CosIntegral}\left (d \left (\frac {e}{f}+x\right )\right ) \sin \left (c-\frac {d e}{f}\right )+f \sin (2 (c+d x))-2 d e \cos \left (c-\frac {d e}{f}\right ) \text {Si}\left (d \left (\frac {e}{f}+x\right )\right )-2 d f x \cos \left (c-\frac {d e}{f}\right ) \text {Si}\left (d \left (\frac {e}{f}+x\right )\right )+2 d e \sin \left (2 c-\frac {2 d e}{f}\right ) \text {Si}\left (\frac {2 d (e+f x)}{f}\right )+2 d f x \sin \left (2 c-\frac {2 d e}{f}\right ) \text {Si}\left (\frac {2 d (e+f x)}{f}\right )}{2 a f^2 (e+f x)} \]

[In]

Integrate[Cos[c + d*x]^3/((e + f*x)^2*(a + a*Sin[c + d*x])),x]

[Out]

(-2*f*Cos[c + d*x] - 2*d*(e + f*x)*Cos[2*c - (2*d*e)/f]*CosIntegral[(2*d*(e + f*x))/f] - 2*d*(e + f*x)*CosInte
gral[d*(e/f + x)]*Sin[c - (d*e)/f] + f*Sin[2*(c + d*x)] - 2*d*e*Cos[c - (d*e)/f]*SinIntegral[d*(e/f + x)] - 2*
d*f*x*Cos[c - (d*e)/f]*SinIntegral[d*(e/f + x)] + 2*d*e*Sin[2*c - (2*d*e)/f]*SinIntegral[(2*d*(e + f*x))/f] +
2*d*f*x*Sin[2*c - (2*d*e)/f]*SinIntegral[(2*d*(e + f*x))/f])/(2*a*f^2*(e + f*x))

Maple [C] (verified)

Result contains complex when optimal does not.

Time = 1.03 (sec) , antiderivative size = 297, normalized size of antiderivative = 1.70

method result size
risch \(\frac {i d \,{\mathrm e}^{-\frac {i \left (c f -d e \right )}{f}} \operatorname {Ei}_{1}\left (i d x +i c -\frac {i \left (c f -d e \right )}{f}\right )}{2 a \,f^{2}}-\frac {i d \,{\mathrm e}^{\frac {i \left (c f -d e \right )}{f}} \operatorname {Ei}_{1}\left (-i d x -i c -\frac {-i c f +i d e}{f}\right )}{2 a \,f^{2}}+\frac {d \,{\mathrm e}^{\frac {2 i \left (c f -d e \right )}{f}} \operatorname {Ei}_{1}\left (-2 i d x -2 i c -\frac {2 \left (-i c f +i d e \right )}{f}\right )}{2 a \,f^{2}}+\frac {d \,{\mathrm e}^{-\frac {2 i \left (c f -d e \right )}{f}} \operatorname {Ei}_{1}\left (2 i d x +2 i c -\frac {2 i \left (c f -d e \right )}{f}\right )}{2 a \,f^{2}}-\frac {\left (-2 d x f -2 d e \right ) \cos \left (d x +c \right )}{2 f a \left (f x +e \right ) \left (-d x f -d e \right )}+\frac {\left (-2 d x f -2 d e \right ) \sin \left (2 d x +2 c \right )}{4 f a \left (f x +e \right ) \left (-d x f -d e \right )}\) \(297\)

[In]

int(cos(d*x+c)^3/(f*x+e)^2/(a+a*sin(d*x+c)),x,method=_RETURNVERBOSE)

[Out]

1/2*I/a*d/f^2*exp(-I*(c*f-d*e)/f)*Ei(1,I*d*x+I*c-I*(c*f-d*e)/f)-1/2*I*d/a/f^2*exp(I*(c*f-d*e)/f)*Ei(1,-I*d*x-I
*c-(-I*c*f+I*d*e)/f)+1/2*d/a/f^2*exp(2*I*(c*f-d*e)/f)*Ei(1,-2*I*d*x-2*I*c-2*(-I*c*f+I*d*e)/f)+1/2/a*d/f^2*exp(
-2*I*(c*f-d*e)/f)*Ei(1,2*I*d*x+2*I*c-2*I*(c*f-d*e)/f)-1/2/f*(-2*d*f*x-2*d*e)/a/(f*x+e)/(-d*f*x-d*e)*cos(d*x+c)
+1/4/f*(-2*d*f*x-2*d*e)/a/(f*x+e)/(-d*f*x-d*e)*sin(2*d*x+2*c)

Fricas [A] (verification not implemented)

none

Time = 0.27 (sec) , antiderivative size = 188, normalized size of antiderivative = 1.07 \[ \int \frac {\cos ^3(c+d x)}{(e+f x)^2 (a+a \sin (c+d x))} \, dx=-\frac {{\left (d f x + d e\right )} \cos \left (-\frac {2 \, {\left (d e - c f\right )}}{f}\right ) \operatorname {Ci}\left (\frac {2 \, {\left (d f x + d e\right )}}{f}\right ) - f \cos \left (d x + c\right ) \sin \left (d x + c\right ) + {\left (d f x + d e\right )} \operatorname {Ci}\left (\frac {d f x + d e}{f}\right ) \sin \left (-\frac {d e - c f}{f}\right ) - {\left (d f x + d e\right )} \sin \left (-\frac {2 \, {\left (d e - c f\right )}}{f}\right ) \operatorname {Si}\left (\frac {2 \, {\left (d f x + d e\right )}}{f}\right ) + {\left (d f x + d e\right )} \cos \left (-\frac {d e - c f}{f}\right ) \operatorname {Si}\left (\frac {d f x + d e}{f}\right ) + f \cos \left (d x + c\right )}{a f^{3} x + a e f^{2}} \]

[In]

integrate(cos(d*x+c)^3/(f*x+e)^2/(a+a*sin(d*x+c)),x, algorithm="fricas")

[Out]

-((d*f*x + d*e)*cos(-2*(d*e - c*f)/f)*cos_integral(2*(d*f*x + d*e)/f) - f*cos(d*x + c)*sin(d*x + c) + (d*f*x +
 d*e)*cos_integral((d*f*x + d*e)/f)*sin(-(d*e - c*f)/f) - (d*f*x + d*e)*sin(-2*(d*e - c*f)/f)*sin_integral(2*(
d*f*x + d*e)/f) + (d*f*x + d*e)*cos(-(d*e - c*f)/f)*sin_integral((d*f*x + d*e)/f) + f*cos(d*x + c))/(a*f^3*x +
 a*e*f^2)

Sympy [F(-1)]

Timed out. \[ \int \frac {\cos ^3(c+d x)}{(e+f x)^2 (a+a \sin (c+d x))} \, dx=\text {Timed out} \]

[In]

integrate(cos(d*x+c)**3/(f*x+e)**2/(a+a*sin(d*x+c)),x)

[Out]

Timed out

Maxima [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 0.33 (sec) , antiderivative size = 309, normalized size of antiderivative = 1.77 \[ \int \frac {\cos ^3(c+d x)}{(e+f x)^2 (a+a \sin (c+d x))} \, dx=-\frac {2 \, d^{2} {\left (E_{2}\left (\frac {i \, d e + i \, {\left (d x + c\right )} f - i \, c f}{f}\right ) + E_{2}\left (-\frac {i \, d e + i \, {\left (d x + c\right )} f - i \, c f}{f}\right )\right )} \cos \left (-\frac {d e - c f}{f}\right ) - d^{2} {\left (-i \, E_{2}\left (\frac {2 \, {\left (-i \, d e - i \, {\left (d x + c\right )} f + i \, c f\right )}}{f}\right ) + i \, E_{2}\left (-\frac {2 \, {\left (-i \, d e - i \, {\left (d x + c\right )} f + i \, c f\right )}}{f}\right )\right )} \cos \left (-\frac {2 \, {\left (d e - c f\right )}}{f}\right ) + 2 \, d^{2} {\left (-i \, E_{2}\left (\frac {i \, d e + i \, {\left (d x + c\right )} f - i \, c f}{f}\right ) + i \, E_{2}\left (-\frac {i \, d e + i \, {\left (d x + c\right )} f - i \, c f}{f}\right )\right )} \sin \left (-\frac {d e - c f}{f}\right ) - d^{2} {\left (E_{2}\left (\frac {2 \, {\left (-i \, d e - i \, {\left (d x + c\right )} f + i \, c f\right )}}{f}\right ) + E_{2}\left (-\frac {2 \, {\left (-i \, d e - i \, {\left (d x + c\right )} f + i \, c f\right )}}{f}\right )\right )} \sin \left (-\frac {2 \, {\left (d e - c f\right )}}{f}\right )}{4 \, {\left (a d e f + {\left (d x + c\right )} a f^{2} - a c f^{2}\right )} d} \]

[In]

integrate(cos(d*x+c)^3/(f*x+e)^2/(a+a*sin(d*x+c)),x, algorithm="maxima")

[Out]

-1/4*(2*d^2*(exp_integral_e(2, (I*d*e + I*(d*x + c)*f - I*c*f)/f) + exp_integral_e(2, -(I*d*e + I*(d*x + c)*f
- I*c*f)/f))*cos(-(d*e - c*f)/f) - d^2*(-I*exp_integral_e(2, 2*(-I*d*e - I*(d*x + c)*f + I*c*f)/f) + I*exp_int
egral_e(2, -2*(-I*d*e - I*(d*x + c)*f + I*c*f)/f))*cos(-2*(d*e - c*f)/f) + 2*d^2*(-I*exp_integral_e(2, (I*d*e
+ I*(d*x + c)*f - I*c*f)/f) + I*exp_integral_e(2, -(I*d*e + I*(d*x + c)*f - I*c*f)/f))*sin(-(d*e - c*f)/f) - d
^2*(exp_integral_e(2, 2*(-I*d*e - I*(d*x + c)*f + I*c*f)/f) + exp_integral_e(2, -2*(-I*d*e - I*(d*x + c)*f + I
*c*f)/f))*sin(-2*(d*e - c*f)/f))/((a*d*e*f + (d*x + c)*a*f^2 - a*c*f^2)*d)

Giac [C] (verification not implemented)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 2.50 (sec) , antiderivative size = 46878, normalized size of antiderivative = 267.87 \[ \int \frac {\cos ^3(c+d x)}{(e+f x)^2 (a+a \sin (c+d x))} \, dx=\text {Too large to display} \]

[In]

integrate(cos(d*x+c)^3/(f*x+e)^2/(a+a*sin(d*x+c)),x, algorithm="giac")

[Out]

-1/2*(d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d
*e/f)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(
1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f
)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*
tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*t
an(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*
c)^4*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan
(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*
x)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*t
an(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*ta
n(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d
*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*
e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*
x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integ
ral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 8*d*f*x*sin_integral
(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_i
ntegral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*imag_part(cos
_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*real_part(
cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*real
_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2
*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*
x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + d*f*x*imag_part(c
os_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - d*f*x*real_part(cos_integral(
2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - d*f*x*real_part(cos_integral(-2*d*x -
2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)
^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*
d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan
(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/
2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(
1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*t
an(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f) + d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2
*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2
*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)
^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*t
an(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*ta
n(1/2*d*e/f)^2 - 8*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d
*e/f)*tan(1/2*d*e/f)^2 - 8*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c
)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 + 2*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan
(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x
)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2
*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 + 6*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/
2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*
x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f)
)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*e*imag_part(cos_integral(-2*d*x -
 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral
(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_inte
gral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 8*d*e*sin_integral(
2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_
integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*imag_part(cos_integral(-d
*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(2*d*x + 2*d*
e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*
tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1
/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4
*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d
*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/
f)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)
^2*tan(1/2*d*e/f)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d
*e/f)^2 + 2*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f) -
 2*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f) + 4*d*f*x
*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f) - 4*d*f*x*imag_part(cos_int
egral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(-
2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 + 2*d*f*x*real_part(cos_integral(d*x + d
*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan
(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*
d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*
c)^4*tan(d*e/f)^2 + d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f
)^2 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - d*e*r
eal_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - 2*d*e*sin_integ
ral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + 2*d*f*x*real_part(cos_integral(d*x
+ d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f) + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f)
)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f) + 4*d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^
2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*
x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*
tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x
)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/
2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*
tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*
e/f)^2*tan(1/2*d*e/f) + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^
3*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*
tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2
*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f
)^2 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 + d*e*im
ag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 - d*e*imag_part(cos
_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(
2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(-2*d*x
- 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(
d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 - 12*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*
x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 12*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/
f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 24*d*f*x*sin_integral(2*(d*f*x + d*e)
/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*e*real_part(cos_integral(2*d*x +
2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*e*real_part(cos_integral(-2
*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 + 2*d*f*x*imag_part(cos_in
tegral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*f*x*imag_part(cos_integral(
-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/
f)*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 + 2*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(
1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/
2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan
(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*
x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*t
an(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)
^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan
(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(
d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*e*real_part(cos_integral(2*d*x + 2*d*e/f)
)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*e*real_part(cos_integral(-2*d*x -
 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_integr
al(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(-2
*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x +
d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*
tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan
(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(
1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(
1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c
)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*
tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)
^2*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*
e/f)^2 - d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e
*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*real_par
t(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*d*e*sin_integral((
d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_integral(d*x + d*e/f
))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1
/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*
x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2
*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^4*
tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*f*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f
*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4 + d*f*x*imag_part(cos_integral(
-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4 + d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x
)^2*tan(1/2*d*x)^2*tan(1/2*c)^4 + d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*ta
n(1/2*c)^4 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4 - 8*d*f*x*real_part(
cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f) - 8*d*f*x*real_part(cos_integ
ral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f) + 2*d*e*imag_part(cos_integral(2*d*x
+ 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f) - 2*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f)
)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f) + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2
*d*x)^2*tan(1/2*c)^4*tan(d*e/f) + 6*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*t
an(1/2*c)^2*tan(d*e/f)^2 + 6*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2
*c)^2*tan(d*e/f)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan
(d*e/f)^2 + 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^
2 + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 + 2*d*e*rea
l_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - 8*d*e*sin_integral(2*
(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - d*f*x*imag_part(cos_integral(d*x + d*e/
f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^4
*tan(d*e/f)^2 - d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - d*f*x*re
al_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - 2*d*f*x*sin_integral((d*f*x + d
*e)/f)*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/
2*c)^4*tan(d*e/f)^2 + d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - d
*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - d*f*x*real_part(cos_i
ntegral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*ta
n(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^
2*tan(1/2*c)^3*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/
2*c)^3*tan(1/2*d*e/f) + 8*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d
*e/f) + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f) + 2*d
*e*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f) + 4*d*f*x*imag_
part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*f*x*ima
g_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*f*x*
sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) + 4*d*f*x*imag_
part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_in
tegral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*f*x*sin_integral((d*f*x + d*e)
/f)*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 4*d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2
*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^
2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^3
*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2
*tan(1/2*d*e/f) + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e
/f) + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d
*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f) - 6*d*f*x*rea
l_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 - 6*d*f*x*real_p
art(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + 4*d*e*imag_part(
cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 4*d*e*imag_part(cos_i
ntegral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integ
ral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(-d*x
- d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 + 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(
d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*ta
n(1/2*c)^4*tan(1/2*d*e/f)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/
f)^2 + d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*f*x*real_pa
rt(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*d*f*x*sin_integral((d*f*x + d*
e)/f)*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan
(1/2*c)^4*tan(1/2*d*e/f)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d
*e/f)^2 + d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*f*x*
real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*d*f*x*sin_integral(
(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 8*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f
))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*f*x*real_part(cos_integral(-2*d*x -
2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 - 12*d*e*imag_part(cos_integral(2*d
*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 12*d*e*imag_part(cos_integ
ral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 24*d*e*sin_integra
l(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*f*x*real_part(co
s_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*f*x*real_part(cos_integ
ral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*f*x*real_part(cos_integral(2*
d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*f*x*real_part(cos_integral(-2*d*
x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 + 2*d*e*imag_part(cos_integral(2*d*x + 2
*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))
*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/
2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 + 2*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^
4*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*ta
n(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)*tan(1
/2*d*e/f)^2 - d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)
^2 + d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f
*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*re
al_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*sin_
integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*e*imag_part(cos_integr
al(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*e*imag_part(cos_
integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_pa
rt(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_
part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 8*d*e*si
n_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*f*x*rea
l_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*f*x*real_par
t(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*f*x*real_part(co
s_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*f*x*real_part(cos
_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*e*imag_part(cos_i
ntegral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*e*imag_part(cos_integral
(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(d*x +
 d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(-d*x - d*e/f))*t
an(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/
2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c
)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^
3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d
*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^
2*tan(1/2*d*e/f)^2 - 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*
e/f)^2 - 4*f*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integr
al(d*x + d*e/f))*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(
1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f
)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f
)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*f*tan(d*x)*tan(1/2*
d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^
2*tan(1/2*d*x)^2*tan(1/2*c)^3 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*ta
n(1/2*c)^3 + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3 + 2*d*f*x*rea
l_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3 + 8*d*f*x*sin_integral(2*(d*f*x + d*
e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x
)^2*tan(1/2*c)^4 + d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4 + d*e*real
_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4 + d*e*real_part(cos_integral(-2*d*
x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4 - 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*
d*x)^2*tan(1/2*c)^4 - 12*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2
*tan(d*e/f) + 12*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*
e/f) - 24*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f) - 8*d*e*real
_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f) - 8*d*e*real_part(cos_i
ntegral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f) + 2*d*f*x*imag_part(cos_integral(
2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f) - 2*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d
*x)^2*tan(1/2*c)^4*tan(d*e/f) + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f) + 2
*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f) - 2*d*f*x*imag_part(cos
_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f) + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*
tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f) + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d
*x)^2*tan(1/2*c)*tan(d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*ta
n(1/2*c)*tan(d*e/f)^2 + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(
d*e/f)^2 + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 8
*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 6*d*e*real_part(cos
_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 + 6*d*e*real_part(cos_integral
(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(2*d*x
 + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)
^2*tan(1/2*c)^3*tan(d*e/f)^2 + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)
^2 + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - 8*d*f*x*sin_integral
(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*ta
n(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(
1/2*c)^3*tan(d*e/f)^2 + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2
+ 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - 8*d*f*x*sin_integra
l(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d
*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)
^2 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - d*e*real_part(cos_int
egral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2
*tan(1/2*c)^4*tan(d*e/f)^2 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2
 + d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - d*e*real_part(cos_inte
gral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - d*e*real_part(cos_integral(-2*d*x - 2*d*e/f)
)*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 - 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^4*t
an(d*e/f)^2 + 2*f*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + 4*d*e*imag_part(cos_integral(d*x + d*e
/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d
*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f) + 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)
^2*tan(1/2*c)^3*tan(1/2*d*e/f) + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*
d*e/f) + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f) + 2*d*f*x*real_p
art(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f) + 2*d*f*x*real_part(cos_integral(-d*
x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f) - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)
^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f) - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1
/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f) + 4*d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*t
an(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*
tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2
*c)*tan(d*e/f)^2*tan(1/2*d*e/f) + 4*d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f
)^2*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*
d*e/f) + 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 4*d*e*imag_
part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_
integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*e*sin_integral((d*f*x +
d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))
*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^4*tan(d*e
/f)^2*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*t
an(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1
/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)
^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 8*d
*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 6*d*e*real_part(c
os_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 - 6*d*e*real_part(cos_in
tegral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos_inte
gral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*
e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1
/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f
)^2 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos
_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(-2*d
*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan
(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1
/2*c)^3*tan(1/2*d*e/f)^2 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^
2 + d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 - d*e*imag_part(cos_inte
gral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*ta
n(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4
*tan(1/2*d*e/f)^2 + 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*e*imag_pa
rt(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 - d*e*imag_part(cos_integral(-d*x -
 d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d
*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^
4*tan(1/2*d*e/f)^2 + 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*f*ta
n(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d
*x)^2*tan(1/2*d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^
2*tan(1/2*d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)
^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(
1/2*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2
*tan(1/2*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 - 12*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2
*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 12*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2
*tan(d*e/f)*tan(1/2*d*e/f)^2 - 24*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan
(1/2*d*e/f)^2 - 12*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1
/2*d*e/f)^2 + 12*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/
2*d*e/f)^2 - 24*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2
- 8*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*e*r
eal_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*e*real_part
(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*e*real_part(cos_
integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 + 2*d*f*x*imag_part(cos_in
tegral(2*d*x + 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*f*x*imag_part(cos_integral(-2*d*x - 2*
d*e/f))*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)^4*tan(d*
e/f)*tan(1/2*d*e/f)^2 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/
2*d*e/f)^2 + d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2
 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*
real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*sin_
integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos_inte
gral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(-2
*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*
e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(
d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c
)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*t
an(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan
(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)
^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan
(1/2*d*e/f)^2 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^
2 + 6*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d
*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*e*rea
l_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*e*real_p
art(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_p
art(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos_integral
(-2*d*x - 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*
tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^3*tan(d*
e/f)^2*tan(1/2*d*e/f)^2 - 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 -
 8*f*tan(d*x)^2*tan(1/2*d*x)*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_integral(d*x + d*e
/f))*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^4*tan(d
*e/f)^2*tan(1/2*d*e/f)^2 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/
f)^2 - d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*d*e*sin_in
tegral((d*f*x + d*e)/f)*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*f*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*
tan(1/2*d*e/f)^2 + 2*f*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 6*d*f*x*real_part(cos_integ
ral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 6*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f
))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*
d*x)^2*tan(1/2*c)^3 - 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3 +
 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3 + 2*d*e*real_part(cos_integ
ral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3 + 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*t
an(1/2*d*x)^2*tan(1/2*c)^3 - d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^4 + d*f*x*imag_p
art(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^4 + d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d
*x)^2*tan(1/2*c)^4 + d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4 - 2*d*f*x*sin_int
egral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^4 - d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan
(1/2*c)^4 + d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4 + d*f*x*real_part(cos_inte
gral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4 + d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d
*x)^2*tan(1/2*c)^4 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^4 + 8*d*f*x*real_part(cos
_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f) + 8*d*f*x*real_part(cos_integral(-
2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f) - 12*d*e*imag_part(cos_integral(2*d*x + 2*d*
e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f) + 12*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan
(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f) - 24*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x
)^2*tan(1/2*c)^2*tan(d*e/f) - 8*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e
/f) - 8*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f) - 8*d*f*x*real_part
(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f) - 8*d*f*x*real_part(cos_integral(-2*d*x
 - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f) + 2*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^
2*tan(1/2*c)^4*tan(d*e/f) - 2*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)
 + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f) + 2*d*e*imag_part(cos_integral(2*d
*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f) - 2*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/
2*d*x)^2*tan(1/2*c)^4*tan(d*e/f) + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f
) + d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2 - d*f*x*imag_part(cos_in
tegral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2 - d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*
tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2 - d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*
x)^2*tan(d*e/f)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2 + 4*d*e*imag_
part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 - 4*d*e*imag_part(cos_in
tegral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 2*d*e*real_part(cos_integral(d*x
 + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan
(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)
^2*tan(1/2*c)*tan(d*e/f)^2 + 6*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/
f)^2 + 6*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 + 6*d*f*x*real_p
art(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 + 6*d*f*x*real_part(cos_integral(-
2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*ta
n(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 + 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*t
an(d*e/f)^2 + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 + 2*d*e*real_par
t(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan
(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3
*tan(d*e/f)^2 + 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 + 2*d
*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 + 2*d*e*real_part(cos_integra
l(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)
^2*tan(1/2*c)^3*tan(d*e/f)^2 - 4*f*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 - d*f*x*imag_part(cos_i
ntegral(d*x + d*e/f))*tan(1/2*c)^4*tan(d*e/f)^2 + d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^4*tan
(d*e/f)^2 - d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f)^2 - d*f*x*real_part(cos_int
egral(-2*d*x - 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f)^2 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*c)^4*tan(d*
e/f)^2 - 2*f*tan(d*x)*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(d*x + d*e/f))*
tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2
*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) + 8*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*ta
n(1/2*c)*tan(1/2*d*e/f) + 4*d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)
- 4*d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f) + 8*d*f*x*sin_integral(
(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f) + 4*d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2
*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^
3*tan(1/2*d*e/f) + 8*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f) + 2*d*e*re
al_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f) + 2*d*e*real_part(cos_integral(-d*x
- d*e/f))*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f) + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*t
an(1/2*c)^4*tan(1/2*d*e/f) + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d
*e/f) - 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f) - 2*d
*e*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f) + 4*d*f*x*imag_
part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_inte
gral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*f*x*sin_integral((d*f*x + d*e)/f)*
tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) + 4*d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^
2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/
2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^
2*tan(1/2*d*e/f) + 4*d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d
*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*f*x*sin_integral((d*
f*x + d*e)/f)*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*c)
^4*tan(d*e/f)^2*tan(1/2*d*e/f) + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2
*d*e/f) - d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + d*f*x*imag_p
art(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + d*f*x*real_part(cos_integral(2*d*
x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan
(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(
1/2*d*e/f)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e
/f)^2 + 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2
- 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*d*e*rea
l_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 8*d*e*sin_integral(
2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 6*d*f*x*real_part(cos_integral(2*d*
x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 - 6*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan
(d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 - 6*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/
2*c)^2*tan(1/2*d*e/f)^2 - 6*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/
2*d*e/f)^2 + 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 4*d*e*i
mag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integr
al(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(d*
x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 + 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e
/f)^2 + 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 4*d*e*im
ag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_int
egral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(-d*x - d*e/f))
*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 + 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c
)^3*tan(1/2*d*e/f)^2 - 4*f*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integ
ral(d*x + d*e/f))*tan(1/2*c)^4*tan(1/2*d*e/f)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^4*tan
(1/2*d*e/f)^2 + d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*f*x*real_part
(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2
*c)^4*tan(1/2*d*e/f)^2 - 2*f*tan(d*x)*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*d*e*imag_part(cos_integ
ral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*e*imag_part(cos_integral(-2*
d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*
tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d
*x)^2*tan(1/2*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*ta
n(1/2*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2
*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*
tan(d*e/f)*tan(1/2*d*e/f)^2 - 12*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/
f)*tan(1/2*d*e/f)^2 + 12*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(
1/2*d*e/f)^2 - 24*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 12
*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 12*d*e
*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 24*d*e*si
n_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*f*x*real_part(cos_
integral(2*d*x + 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*f*x*real_part(cos_integral(-2*d*x -
2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 + 2*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)
^4*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f)*tan(1
/2*d*e/f)^2 + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)^4*tan(d*e/f)*tan(1/2*d*e/f)^2 - d*f*x*imag_part
(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integral(-d*x - d*e
/f))*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(
d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*
e/f)^2 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*imag_part(cos_
integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integral(-d*x - d*e/
f))*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)
^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2
*tan(1/2*d*e/f)^2 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*e
*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*e*imag_par
t(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_in
tegral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(-d*x -
 d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^
2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan
(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2
*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d
*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*
tan(1/2*d*e/f)^2 + 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)
^2 + 4*f*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*f*x*real_part(cos_integral(2
*d*x + 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f)
)*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 12*f*tan(d*x)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*
d*e/f)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*e*i
mag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_inte
gral(d*x + d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(-d*x - d*e/f))*ta
n(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)^3*tan(d*e/f)^2*tan
(1/2*d*e/f)^2 - 4*f*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*f*tan(1/2*d*x)^2*tan(1/2*c)^3*ta
n(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*f*tan(d*x)*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_
integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c) + 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*
e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c) + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d
*x)^2*tan(1/2*c) + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c) - 8*d*f*
x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c) - 6*d*e*real_part(cos_integral(2*d*x +
2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 6*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*
tan(1/2*d*x)^2*tan(1/2*c)^2 + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3 - 4*d*f
*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3 + 2*d*f*x*real_part(cos_integral(d*x + d*
e/f))*tan(d*x)^2*tan(1/2*c)^3 + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^3 + 8*d*f*
x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^3 + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*t
an(1/2*d*x)^2*tan(1/2*c)^3 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3 + 2
*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3 + 2*d*f*x*real_part(cos_integral(-d*x
- d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^3 -
d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)^4 + d*e*imag_part(cos_integral(-d*x - d*e/f))*t
an(d*x)^2*tan(1/2*c)^4 + d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4 + d*e*real_part(
cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^4 - 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1
/2*c)^4 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4 + d*e*imag_part(cos_integral(-d
*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4 + d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2
*c)^4 + d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^4 - 2*d*e*sin_integral((d*f*x
+ d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^4 + 2*f*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4 + 2*d*f*x*imag_part(cos_int
egral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f) - 2*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f
))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f) + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*ta
n(d*e/f) + 8*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f) + 8*
d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f) - 12*d*f*x*imag_
part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f) + 12*d*f*x*imag_part(cos_integral(-2*d*
x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f) - 24*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2
*c)^2*tan(d*e/f) - 12*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f) +
12*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f) - 24*d*f*x*sin_integ
ral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f) - 8*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))
*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f) - 8*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3*
tan(d*e/f) - 8*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f) - 8*d*e*rea
l_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(d*e/f) + 2*d*f*x*imag_part(cos_integral
(2*d*x + 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f) - 2*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^4*ta
n(d*e/f) + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)^4*tan(d*e/f) + d*e*imag_part(cos_integral(d*x +
d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2 - d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*
d*x)^2*tan(d*e/f)^2 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2 - d*
e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2 + 2*d*e*sin_integral((d*f*x
 + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^
2*tan(1/2*c)*tan(d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)
^2 + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 2*d*f*x*real_part(cos_i
ntegral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2
*tan(1/2*c)*tan(d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/
f)^2 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 2*d*f*x*real
_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 2*d*f*x*real_part(cos_integral(-d*x
- d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(
1/2*c)*tan(d*e/f)^2 + 6*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 + 6*
d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 + 6*d*e*real_part(cos_integ
ral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 + 6*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f
))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^3*ta
n(d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2 + 2*d*f*x*real_part(c
os_integral(d*x + d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2 + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)
^3*tan(d*e/f)^2 - 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)^3*tan(d*e/f)^2 - 8*f*tan(d*x)^2*tan(1/2*d
*x)*tan(1/2*c)^3*tan(d*e/f)^2 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*c)^4*tan(d*e/f)^2 + d*e*imag_
part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^4*tan(d*e/f)^2 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(
1/2*c)^4*tan(d*e/f)^2 - d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f)^2 - 2*d*e*sin_in
tegral((d*f*x + d*e)/f)*tan(1/2*c)^4*tan(d*e/f)^2 - 2*f*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + 2*f*tan(1/2*d*x
)^2*tan(1/2*c)^4*tan(d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2
*d*e/f) - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f) + 4*d*e*imag_
part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_inte
gral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) + 8*d*e*sin_integral((d*f*x + d*e)/f)*
tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) + 4*d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*ta
n(1/2*c)^3*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)
 + 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f) + 4*d*e*imag_part(cos_integral(d
*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*
d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f) + 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^3*tan(1/2*d
*e/f) + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(1/2*c)^4*tan(1/2*d*e/f) + 2*d*f*x*real_part(cos_integ
ral(-d*x - d*e/f))*tan(1/2*c)^4*tan(1/2*d*e/f) - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(d
*e/f)^2*tan(1/2*d*e/f) - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)
- 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f) - 2*d*f*x*real_part(
cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f) + 4*d*e*imag_part(cos_integral(d*x + d*
e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)
^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)*tan(d*e/
f)^2*tan(1/2*d*e/f) + 4*d*e*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/
2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) +
 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) + 4*d*e*imag_part(c
os_integral(d*x + d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f
))*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*c)^3*tan(d*e/f)^2*ta
n(1/2*d*e/f) - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + d*e*imag_
part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(2*d*x
 + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*
x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 - 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d
*e/f)^2 - 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 + 4*d*f*x*im
ag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integra
l(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x
)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)
^2 - 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 + 4*d*f*x*ima
g_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_inte
gral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*
tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)
*tan(1/2*d*e/f)^2 - 6*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 -
6*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 - 6*d*e*real_part(cos
_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 - 6*d*e*real_part(cos_integral(-2*d*x
 - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*t
an(1/2*c)^3*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^3*tan(1/2*d*e/f)^2
 - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral
(-d*x - d*e/f))*tan(1/2*c)^3*tan(1/2*d*e/f)^2 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)^3*tan(1/2*d
*e/f)^2 - 8*f*tan(d*x)^2*tan(1/2*d*x)*tan(1/2*c)^3*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_integral(d*x + d*e/f))
*tan(1/2*c)^4*tan(1/2*d*e/f)^2 - d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*e
*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(-2*d*x -
2*d*e/f))*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*c)^4*tan(1/2*d*e/f)^2 -
2*f*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*f*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*d*f*x*imag
_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*f*x*imag_part(cos_integral(-
2*d*x - 2*d*e/f))*tan(d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*
tan(d*e/f)*tan(1/2*d*e/f)^2 + 2*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)*tan(1
/2*d*e/f)^2 - 2*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4
*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*e*real_part(cos_integr
al(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*e*real_part(cos_integral(-2*d*x -
 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*
tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/
2*d*x)^2*tan(1/2*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 - 12*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)
^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 12*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f)*ta
n(1/2*d*e/f)^2 - 24*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*e*rea
l_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 - 8*d*e*real_part(cos_integral(
-2*d*x - 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)*tan(1/2*d*e/f)^2 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x
)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*
d*e/f)^2 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*real_pa
rt(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*sin_integral((d*f*x + d*e)
/f)*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(d*e
/f)^2*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^
2 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*real_part(
cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*sin_integral((d*f*x + d*e
)/f)*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*f*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^
2 + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_p
art(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(
d*x + d*e/f))*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2
*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*
e/f)^2 - 8*f*tan(d*x)^2*tan(1/2*d*x)*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*e*real_part(cos_integral(2
*d*x + 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 6*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*
tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 8*f*tan(1/2*d*x)*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*f
*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*
x)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2 + d*f*x*real_part(cos_integral(2*
d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2 + d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2
*d*x)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2 - 4*d*e*imag_part(cos_integral(2*d*x
 + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c) + 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2
*tan(1/2*d*x)^2*tan(1/2*c) + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c) +
 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c) - 8*d*e*sin_integral(2*(d*f*
x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c) - 6*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2
*tan(1/2*c)^2 - 6*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2 - 6*d*f*x*real_part(
cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 6*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))
*tan(1/2*d*x)^2*tan(1/2*c)^2 + 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3 - 4*d*e*
imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^3 + 2*d*e*real_part(cos_integral(d*x + d*e/f))
*tan(d*x)^2*tan(1/2*c)^3 + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)^3 + 8*d*e*sin_int
egral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^3 + 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)
^2*tan(1/2*c)^3 - 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3 + 2*d*e*real_par
t(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^3 + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2
*d*x)^2*tan(1/2*c)^3 + 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)^3 - 4*f*tan(d*x)^2*tan(
1/2*d*x)^2*tan(1/2*c)^3 - d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*c)^4 + d*f*x*imag_part(cos_integr
al(-d*x - d*e/f))*tan(1/2*c)^4 + d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^4 + d*f*x*real_part
(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^4 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*c)^4 - 2*f*tan(d
*x)*tan(1/2*d*x)^2*tan(1/2*c)^4 + 2*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan
(d*e/f) - 2*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f) + 4*d*e*sin_int
egral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f) + 8*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f
))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f) + 8*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*
tan(d*e/f) + 8*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f) + 8*d*f*x*r
eal_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f) - 12*d*e*imag_part(cos_integral(
2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f) + 12*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*
x)^2*tan(1/2*c)^2*tan(d*e/f) - 24*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f) - 12*
d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f) + 12*d*e*imag_part(cos_int
egral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f) - 24*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1
/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f) - 8*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f) -
 8*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f) + 2*d*e*imag_part(cos_integral(2*d*
x + 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f) - 2*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^4*tan(d*e/f
) + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)^4*tan(d*e/f) + d*f*x*imag_part(cos_integral(d*x + d*e/f))
*tan(d*x)^2*tan(d*e/f)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(d*e/f)^2 - d*f*x*real_pa
rt(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(d*e/f)^2 - d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*ta
n(d*x)^2*tan(d*e/f)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(d*e/f)^2 + d*f*x*imag_part(cos_in
tegral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*
tan(d*e/f)^2 - d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2 - d*f*x*real_part(co
s_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)
^2*tan(d*e/f)^2 + 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2 - 4*d*e*im
ag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 2*d*e*real_part(cos_integral(d*x
+ d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*
c)*tan(d*e/f)^2 + 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 4*d*e*imag_part(c
os_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 - 4*d*e*imag_part(cos_integral(-2*d*x - 2
*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*ta
n(1/2*c)*tan(d*e/f)^2 + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 8
*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 4*f*tan(d*x)^2*tan(1/2*d*x)^2*ta
n(1/2*c)*tan(d*e/f)^2 + 6*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f)^2 + 6*d*f*x*r
eal_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f)^2 + 12*f*tan(d*x)*tan(1/2*d*x)^2*tan(1/2*c)^2
*tan(d*e/f)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2 + 4*d*e*imag_part(cos
_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2 + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*c)
^3*tan(d*e/f)^2 + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^3*tan(d*e/f)^2 - 8*d*e*sin_integral(2
*(d*f*x + d*e)/f)*tan(1/2*c)^3*tan(d*e/f)^2 - 4*f*tan(d*x)^2*tan(1/2*c)^3*tan(d*e/f)^2 + 4*f*tan(1/2*d*x)^2*ta
n(1/2*c)^3*tan(d*e/f)^2 - 2*f*tan(d*x)*tan(1/2*c)^4*tan(d*e/f)^2 - 2*d*e*real_part(cos_integral(d*x + d*e/f))*
tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f) - 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)
^2*tan(1/2*d*e/f) + 4*d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) - 4*d*f*
x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) + 8*d*f*x*sin_integral((d*f*x + d
*e)/f)*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) + 4*d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(
1/2*c)*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)
 + 8*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) + 4*d*f*x*imag_part(cos_inte
gral(d*x + d*e/f))*tan(1/2*c)^3*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^3*ta
n(1/2*d*e/f) + 8*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*c)^3*tan(1/2*d*e/f) + 2*d*e*real_part(cos_integra
l(d*x + d*e/f))*tan(1/2*c)^4*tan(1/2*d*e/f) + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^4*tan(1/2
*d*e/f) - 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f) - 2*d*e*real_part(
cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f) - 2*d*e*real_part(cos_integral(d*x + d*e/f)
)*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f) - 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(
d*e/f)^2*tan(1/2*d*e/f) + 4*d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)
- 4*d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*f*x*sin_integral(
(d*f*x + d*e)/f)*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) - d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^
2*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*e/f)^2 + d*f*x*real_part
(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*e/f)^2 + d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*
tan(d*x)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*e/f)^2 - d*f*x*imag_p
art(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integral(-d*x - d*e/f))*t
an(1/2*d*x)^2*tan(1/2*d*e/f)^2 + d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*d*e/f)^
2 + d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*sin_integral((d*
f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(
1/2*c)*tan(1/2*d*e/f)^2 + 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)
^2 - 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_i
ntegral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)
^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(
1/2*d*e/f)^2 + 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*
d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_inte
gral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2
*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 + 4*f*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 6*d*f*x*real
_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^2*tan(1/2*d*e/f)^2 - 6*d*f*x*real_part(cos_integral(-2*d*x - 2
*d*e/f))*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + 12*f*tan(d*x)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + 4*d*e*im
ag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 4*d*e*imag_part(cos_integral(-2*d*x - 2
*d*e/f))*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*c)^3*tan(1/2*d*e/f
)^2 - 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^3*tan(1/2*d*e/f)^2 + 8*d*e*sin_integral(2*(d*f*x
+ d*e)/f)*tan(1/2*c)^3*tan(1/2*d*e/f)^2 - 4*f*tan(d*x)^2*tan(1/2*c)^3*tan(1/2*d*e/f)^2 + 4*f*tan(1/2*d*x)^2*ta
n(1/2*c)^3*tan(1/2*d*e/f)^2 - 2*f*tan(d*x)*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*d*e*imag_part(cos_integral(2*d*x
+ 2*d*e/f))*tan(d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^
2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 +
 2*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*e*imag_part(c
os_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*e*sin_integral(2*(d*f*x + d*e)
/f)*tan(1/2*d*x)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)*t
an(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)*tan(d*e/f)*tan(1/2*d
*e/f)^2 - 12*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 + 12*d*e*im
ag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - 24*d*e*sin_integral(2*(d*f*
x + d*e)/f)*tan(1/2*c)^2*tan(d*e/f)*tan(1/2*d*e/f)^2 - d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*e/f)^2
*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*real_par
t(cos_integral(2*d*x + 2*d*e/f))*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f
))*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*f*t
an(d*x)*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*
c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)*tan(d*e/f)^2*tan
(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*re
al_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 8*d*e*sin_integral(2*(d*f*x + d
*e)/f)*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 4*f*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*
f*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 12*f*tan(d*x)*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*
e/f)^2 + 4*f*tan(1/2*c)^3*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*
tan(1/2*d*x)^2 - d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2 + d*e*real_part(cos_integ
ral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*x)^2 + d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan
(1/2*d*x)^2 + 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*d*x)^2 - 4*d*f*x*imag_part(cos_integral(2
*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c) + 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*
c) + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c) + 2*d*f*x*real_part(cos_integral(-d*x
- d*e/f))*tan(d*x)^2*tan(1/2*c) - 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c) - 4*d*f*x*imag
_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c) + 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e
/f))*tan(1/2*d*x)^2*tan(1/2*c) + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c) + 2*d*
f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c) - 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*
tan(1/2*d*x)^2*tan(1/2*c) - 6*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2 - 6*d*e*rea
l_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)^2 - 6*d*e*real_part(cos_integral(2*d*x + 2*d*e/f)
)*tan(1/2*d*x)^2*tan(1/2*c)^2 - 6*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)^2 +
4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^3 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/
f))*tan(1/2*c)^3 + 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(1/2*c)^3 + 2*d*f*x*real_part(cos_integral(
-d*x - d*e/f))*tan(1/2*c)^3 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)^3 - 8*f*tan(d*x)^2*tan(1/2*d*
x)*tan(1/2*c)^3 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*c)^4 + d*e*imag_part(cos_integral(-d*x - d*
e/f))*tan(1/2*c)^4 + d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^4 + d*e*real_part(cos_integral(-2
*d*x - 2*d*e/f))*tan(1/2*c)^4 - 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*c)^4 - 2*f*tan(d*x)^2*tan(1/2*c)^4
 + 2*f*tan(1/2*d*x)^2*tan(1/2*c)^4 + 2*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(d*e/f) -
2*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(d*e/f) + 4*d*f*x*sin_integral(2*(d*f*x + d*e)
/f)*tan(d*x)^2*tan(d*e/f) + 2*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f) - 2*d*f
*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f) + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/
f)*tan(1/2*d*x)^2*tan(d*e/f) + 8*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)
 + 8*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(d*e/f) + 8*d*e*real_part(cos_inte
gral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f) + 8*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*
tan(1/2*d*x)^2*tan(1/2*c)*tan(d*e/f) - 12*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^2*tan(d*e/
f) + 12*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f) - 24*d*f*x*sin_integral(2*(d*f
*x + d*e)/f)*tan(1/2*c)^2*tan(d*e/f) - 8*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f)
- 8*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^3*tan(d*e/f) + d*e*imag_part(cos_integral(d*x + d
*e/f))*tan(d*x)^2*tan(d*e/f)^2 - d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(d*e/f)^2 - d*e*real_
part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(d*e/f)^2 - d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*ta
n(d*x)^2*tan(d*e/f)^2 + 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan(d*e/f)^2 + d*e*imag_part(cos_integr
al(d*x + d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2 - d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(d*
e/f)^2 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2 - d*e*real_part(cos_integral
(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f)^2 + 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(d*e/f
)^2 - 2*f*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2 + 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)
*tan(d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)*tan(d*e/f)^2 + 2*d*f*x*real_part(
cos_integral(d*x + d*e/f))*tan(1/2*c)*tan(d*e/f)^2 + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)*
tan(d*e/f)^2 + 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)*tan(d*e/f)^2 - 8*f*tan(d*x)^2*tan(1/2*d*x)*t
an(1/2*c)*tan(d*e/f)^2 + 6*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f)^2 + 6*d*e*real
_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f)^2 - 8*f*tan(1/2*d*x)*tan(1/2*c)^3*tan(d*e/f)^2 -
 2*f*tan(1/2*c)^4*tan(d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*e/f) - 2*d*
f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*d*e/f) - 2*d*f*x*real_part(cos_integral(d*x + d*e
/f))*tan(1/2*d*x)^2*tan(1/2*d*e/f) - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*d*e/
f) + 4*d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_int
egral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) + 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2*tan
(1/2*c)*tan(1/2*d*e/f) + 4*d*e*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) -
 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) + 8*d*e*sin_integral((d*
f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f) + 4*d*e*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*c)
^3*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^3*tan(1/2*d*e/f) + 8*d*e*sin_integr
al((d*f*x + d*e)/f)*tan(1/2*c)^3*tan(1/2*d*e/f) - 2*d*f*x*real_part(cos_integral(d*x + d*e/f))*tan(d*e/f)^2*ta
n(1/2*d*e/f) - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(d*e/f)^2*tan(1/2*d*e/f) + 4*d*e*imag_part(cos
_integral(d*x + d*e/f))*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*t
an(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f) + 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*
e/f) - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_integral(-d*x
- d*e/f))*tan(d*x)^2*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*d*e/f)
^2 + d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*d*e/f)^2 - 2*d*e*sin_integral((d*f*x + d
*e)/f)*tan(d*x)^2*tan(1/2*d*e/f)^2 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2
+ d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(2*d*x
 + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*ta
n(1/2*d*e/f)^2 - 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 - 2*f*tan(d*x)^2*tan(1/2*
d*x)^2*tan(1/2*d*e/f)^2 - 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)*tan(1/2*d*e/f)^2 + 4*d*f
*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x
+ d*e/f))*tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)*tan(1/2*d*e/f
)^2 - 8*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)*tan(1/2*d*e/f)^2 - 8*f*tan(d*x)^2*tan(1/2*d*x)*tan(1/
2*c)*tan(1/2*d*e/f)^2 - 6*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^2*tan(1/2*d*e/f)^2 - 6*d*e*r
eal_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^2*tan(1/2*d*e/f)^2 - 8*f*tan(1/2*d*x)*tan(1/2*c)^3*tan(1/2
*d*e/f)^2 - 2*f*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*e/f)*ta
n(1/2*d*e/f)^2 - 2*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*f*x*sin_i
ntegral(2*(d*f*x + d*e)/f)*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/
2*c)*tan(d*e/f)*tan(1/2*d*e/f)^2 + 8*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)*tan(d*e/f)*tan(1
/2*d*e/f)^2 - d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_integ
ral(-d*x - d*e/f))*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*e/f)^2*t
an(1/2*d*e/f)^2 - d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*d*e*sin_inte
gral((d*f*x + d*e)/f)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*f*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 2*f*tan(1
/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 8*f*tan(1/2*d*x)*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*f*x*im
ag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2 + d*f*x
*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2 + d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x
)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2 + d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*
x)^2 - d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2 + d*f*x*real_part(cos_integral(2*d*x + 2*d*e
/f))*tan(1/2*d*x)^2 + d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2 + 2*d*f*x*sin_integral((d
*f*x + d*e)/f)*tan(1/2*d*x)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(1/2*c) + 4*d*e*i
mag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(1/2*c) + 2*d*e*real_part(cos_integral(d*x + d*e/f))*ta
n(d*x)^2*tan(1/2*c) + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2*tan(1/2*c) - 8*d*e*sin_integral(2
*(d*f*x + d*e)/f)*tan(d*x)^2*tan(1/2*c) - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(1/
2*c) + 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2*tan(1/2*c) + 2*d*e*real_part(cos_integra
l(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*c) + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2
*c) - 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(1/2*c) + 4*f*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*
c) - 6*d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^2 - 6*d*f*x*real_part(cos_integral(-2*d*x - 2
*d*e/f))*tan(1/2*c)^2 + 12*f*tan(d*x)*tan(1/2*d*x)^2*tan(1/2*c)^2 + 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e
/f))*tan(1/2*c)^3 - 4*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^3 + 2*d*e*real_part(cos_integra
l(d*x + d*e/f))*tan(1/2*c)^3 + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)^3 + 8*d*e*sin_integral(2
*(d*f*x + d*e)/f)*tan(1/2*c)^3 - 4*f*tan(d*x)^2*tan(1/2*c)^3 + 4*f*tan(1/2*d*x)^2*tan(1/2*c)^3 - 2*f*tan(d*x)*
tan(1/2*c)^4 + 2*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2*tan(d*e/f) - 2*d*e*imag_part(cos_inte
gral(-2*d*x - 2*d*e/f))*tan(d*x)^2*tan(d*e/f) + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*x)^2*tan(d*e/f) +
2*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f) - 2*d*e*imag_part(cos_integral(-2*d*x
 - 2*d*e/f))*tan(1/2*d*x)^2*tan(d*e/f) + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*d*x)^2*tan(d*e/f) + 8*d
*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)*tan(d*e/f) + 8*d*f*x*real_part(cos_integral(-2*d*x -
2*d*e/f))*tan(1/2*c)*tan(d*e/f) - 12*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f) + 12
*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^2*tan(d*e/f) - 24*d*e*sin_integral(2*(d*f*x + d*e)/f
)*tan(1/2*c)^2*tan(d*e/f) + d*f*x*imag_part(cos_integral(d*x + d*e/f))*tan(d*e/f)^2 - d*f*x*imag_part(cos_inte
gral(-d*x - d*e/f))*tan(d*e/f)^2 - d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*e/f)^2 - d*f*x*real_pa
rt(cos_integral(-2*d*x - 2*d*e/f))*tan(d*e/f)^2 + 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(d*e/f)^2 - 2*f*tan
(d*x)*tan(1/2*d*x)^2*tan(d*e/f)^2 + 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)*tan(d*e/f)^2 - 4
*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)*tan(d*e/f)^2 + 2*d*e*real_part(cos_integral(d*x + d*
e/f))*tan(1/2*c)*tan(d*e/f)^2 + 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)*tan(d*e/f)^2 + 8*d*e*si
n_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)*tan(d*e/f)^2 + 4*f*tan(d*x)^2*tan(1/2*c)*tan(d*e/f)^2 - 4*f*tan(1/2*d
*x)^2*tan(1/2*c)*tan(d*e/f)^2 + 12*f*tan(d*x)*tan(1/2*c)^2*tan(d*e/f)^2 + 4*f*tan(1/2*c)^3*tan(d*e/f)^2 - 2*d*
e*real_part(cos_integral(d*x + d*e/f))*tan(d*x)^2*tan(1/2*d*e/f) - 2*d*e*real_part(cos_integral(-d*x - d*e/f))
*tan(d*x)^2*tan(1/2*d*e/f) - 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x)^2*tan(1/2*d*e/f) - 2*d*e*
real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2*tan(1/2*d*e/f) + 4*d*f*x*imag_part(cos_integral(d*x + d*e
/f))*tan(1/2*c)*tan(1/2*d*e/f) - 4*d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)*tan(1/2*d*e/f) + 8*d
*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*c)*tan(1/2*d*e/f) - 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(
d*e/f)^2*tan(1/2*d*e/f) - 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(d*e/f)^2*tan(1/2*d*e/f) - d*f*x*imag
_part(cos_integral(d*x + d*e/f))*tan(1/2*d*e/f)^2 + d*f*x*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*e/f)
^2 + d*f*x*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*d*e/f)^2 + d*f*x*real_part(cos_integral(-2*d*x - 2
*d*e/f))*tan(1/2*d*e/f)^2 - 2*d*f*x*sin_integral((d*f*x + d*e)/f)*tan(1/2*d*e/f)^2 - 2*f*tan(d*x)*tan(1/2*d*x)
^2*tan(1/2*d*e/f)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c)*tan(1/2*d*e/f)^2 + 4*d*e*imag_
part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(d*x + d*e/f))*
tan(1/2*c)*tan(1/2*d*e/f)^2 - 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)*tan(1/2*d*e/f)^2 - 8*d*e*
sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c)*tan(1/2*d*e/f)^2 + 4*f*tan(d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 - 4*f
*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*d*e/f)^2 + 12*f*tan(d*x)*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + 4*f*tan(1/2*c)^3*t
an(1/2*d*e/f)^2 + 2*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*e/f)*tan(1/2*d*e/f)^2 - 2*d*e*imag_part
(cos_integral(-2*d*x - 2*d*e/f))*tan(d*e/f)*tan(1/2*d*e/f)^2 + 4*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(d*e/f
)*tan(1/2*d*e/f)^2 - 2*f*tan(d*x)*tan(d*e/f)^2*tan(1/2*d*e/f)^2 - 4*f*tan(1/2*c)*tan(d*e/f)^2*tan(1/2*d*e/f)^2
 + d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*x)^2 - d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(d*x)^2
+ d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*x)^2 + d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(
d*x)^2 + 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*x)^2 + d*e*imag_part(cos_integral(d*x + d*e/f))*tan(1/2*d*x
)^2 - d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*x)^2 + d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*
tan(1/2*d*x)^2 + d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*x)^2 + 2*d*e*sin_integral((d*f*x + d*
e)/f)*tan(1/2*d*x)^2 - 2*f*tan(d*x)^2*tan(1/2*d*x)^2 - 4*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/
2*c) + 4*d*f*x*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c) + 2*d*f*x*real_part(cos_integral(d*x + d*e
/f))*tan(1/2*c) + 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*c) - 8*d*f*x*sin_integral(2*(d*f*x + d
*e)/f)*tan(1/2*c) - 8*f*tan(d*x)^2*tan(1/2*d*x)*tan(1/2*c) - 6*d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*ta
n(1/2*c)^2 - 6*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)^2 - 8*f*tan(1/2*d*x)*tan(1/2*c)^3 - 2*
f*tan(1/2*c)^4 + 2*d*f*x*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*e/f) - 2*d*f*x*imag_part(cos_integral(
-2*d*x - 2*d*e/f))*tan(d*e/f) + 4*d*f*x*sin_integral(2*(d*f*x + d*e)/f)*tan(d*e/f) + 8*d*e*real_part(cos_integ
ral(2*d*x + 2*d*e/f))*tan(1/2*c)*tan(d*e/f) + 8*d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*c)*tan(d
*e/f) + d*e*imag_part(cos_integral(d*x + d*e/f))*tan(d*e/f)^2 - d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(
d*e/f)^2 - d*e*real_part(cos_integral(2*d*x + 2*d*e/f))*tan(d*e/f)^2 - d*e*real_part(cos_integral(-2*d*x - 2*d
*e/f))*tan(d*e/f)^2 + 2*d*e*sin_integral((d*f*x + d*e)/f)*tan(d*e/f)^2 + 2*f*tan(d*x)^2*tan(d*e/f)^2 - 2*f*tan
(1/2*d*x)^2*tan(d*e/f)^2 - 8*f*tan(1/2*d*x)*tan(1/2*c)*tan(d*e/f)^2 - 2*d*f*x*real_part(cos_integral(d*x + d*e
/f))*tan(1/2*d*e/f) - 2*d*f*x*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*e/f) + 4*d*e*imag_part(cos_integ
ral(d*x + d*e/f))*tan(1/2*c)*tan(1/2*d*e/f) - 4*d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*c)*tan(1/2*d
*e/f) + 8*d*e*sin_integral((d*f*x + d*e)/f)*tan(1/2*c)*tan(1/2*d*e/f) - d*e*imag_part(cos_integral(d*x + d*e/f
))*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(
2*d*x + 2*d*e/f))*tan(1/2*d*e/f)^2 + d*e*real_part(cos_integral(-2*d*x - 2*d*e/f))*tan(1/2*d*e/f)^2 - 2*d*e*si
n_integral((d*f*x + d*e)/f)*tan(1/2*d*e/f)^2 + 2*f*tan(d*x)^2*tan(1/2*d*e/f)^2 - 2*f*tan(1/2*d*x)^2*tan(1/2*d*
e/f)^2 - 8*f*tan(1/2*d*x)*tan(1/2*c)*tan(1/2*d*e/f)^2 + 2*f*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + d*f*x*imag_part(co
s_integral(d*x + d*e/f)) - d*f*x*imag_part(cos_integral(-d*x - d*e/f)) + d*f*x*real_part(cos_integral(2*d*x +
2*d*e/f)) + d*f*x*real_part(cos_integral(-2*d*x - 2*d*e/f)) + 2*d*f*x*sin_integral((d*f*x + d*e)/f) - 2*f*tan(
d*x)*tan(1/2*d*x)^2 - 4*d*e*imag_part(cos_integral(2*d*x + 2*d*e/f))*tan(1/2*c) + 4*d*e*imag_part(cos_integral
(-2*d*x - 2*d*e/f))*tan(1/2*c) + 2*d*e*real_part(cos_integral(d*x + d*e/f))*tan(1/2*c) + 2*d*e*real_part(cos_i
ntegral(-d*x - d*e/f))*tan(1/2*c) - 8*d*e*sin_integral(2*(d*f*x + d*e)/f)*tan(1/2*c) + 4*f*tan(d*x)^2*tan(1/2*
c) - 4*f*tan(1/2*d*x)^2*tan(1/2*c) + 12*f*tan(d*x)*tan(1/2*c)^2 + 4*f*tan(1/2*c)^3 + 2*d*e*imag_part(cos_integ
ral(2*d*x + 2*d*e/f))*tan(d*e/f) - 2*d*e*imag_part(cos_integral(-2*d*x - 2*d*e/f))*tan(d*e/f) + 4*d*e*sin_inte
gral(2*(d*f*x + d*e)/f)*tan(d*e/f) - 2*f*tan(d*x)*tan(d*e/f)^2 - 4*f*tan(1/2*c)*tan(d*e/f)^2 - 2*d*e*real_part
(cos_integral(d*x + d*e/f))*tan(1/2*d*e/f) - 2*d*e*real_part(cos_integral(-d*x - d*e/f))*tan(1/2*d*e/f) - 2*f*
tan(d*x)*tan(1/2*d*e/f)^2 - 4*f*tan(1/2*c)*tan(1/2*d*e/f)^2 + d*e*imag_part(cos_integral(d*x + d*e/f)) - d*e*i
mag_part(cos_integral(-d*x - d*e/f)) + d*e*real_part(cos_integral(2*d*x + 2*d*e/f)) + d*e*real_part(cos_integr
al(-2*d*x - 2*d*e/f)) + 2*d*e*sin_integral((d*f*x + d*e)/f) + 2*f*tan(d*x)^2 - 2*f*tan(1/2*d*x)^2 - 8*f*tan(1/
2*d*x)*tan(1/2*c) + 2*f*tan(d*e/f)^2 + 2*f*tan(1/2*d*e/f)^2 - 2*f*tan(d*x) - 4*f*tan(1/2*c) + 2*f)/(a*f^3*x*ta
n(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*e*f^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2
*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*f^3*x*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + a*f^3*x*ta
n(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*a*f^3*x*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d
*e/f)^2*tan(1/2*d*e/f)^2 + a*f^3*x*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*f^3*x*tan(1/2*d*x
)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*e*f^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 +
 a*e*f^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + 2*a*e*f^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2
*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*e*f^2*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*e*f^2*
tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*f^3*x*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4 + 2
*a*f^3*x*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 + a*f^3*x*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 +
a*f^3*x*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + 2*a*f^3*x*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*d*
e/f)^2 + a*f^3*x*tan(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + a*f^3*x*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)
^2 + a*f^3*x*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*a*f^3*x*tan(d*x)^2*tan(1/2*c)^2*tan(d
*e/f)^2*tan(1/2*d*e/f)^2 + 2*a*f^3*x*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*f^3*x*tan(1
/2*c)^4*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*e*f^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^4 + 2*a*e*f^2*tan(d*x)^2*
tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 + a*e*f^2*tan(d*x)^2*tan(1/2*c)^4*tan(d*e/f)^2 + a*e*f^2*tan(1/2*d*x)
^2*tan(1/2*c)^4*tan(d*e/f)^2 + 2*a*e*f^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + a*e*f^2*tan
(d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + a*e*f^2*tan(1/2*d*x)^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + a*e*f^2*tan(d*x
)^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*a*e*f^2*tan(d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/
f)^2 + 2*a*e*f^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*e*f^2*tan(1/2*c)^4*tan(d*e/f)^2
*tan(1/2*d*e/f)^2 + 2*a*f^3*x*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2 + a*f^3*x*tan(d*x)^2*tan(1/2*c)^4 + a*f^3
*x*tan(1/2*d*x)^2*tan(1/2*c)^4 + a*f^3*x*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2 + 2*a*f^3*x*tan(d*x)^2*tan(1/2
*c)^2*tan(d*e/f)^2 + 2*a*f^3*x*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 + a*f^3*x*tan(1/2*c)^4*tan(d*e/f)^2 +
a*f^3*x*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + 2*a*f^3*x*tan(d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + 2*a*
f^3*x*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + a*f^3*x*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + a*f^3*x*tan(d*x)^
2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*f^3*x*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*a*f^3*x*tan(1/2*c)^
2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*a*e*f^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2 + a*e*f^2*tan(d*x)^2*tan(1/
2*c)^4 + a*e*f^2*tan(1/2*d*x)^2*tan(1/2*c)^4 + a*e*f^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(d*e/f)^2 + 2*a*e*f^2*tan(
d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 + 2*a*e*f^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(d*e/f)^2 + a*e*f^2*tan(1/2*c)^4*t
an(d*e/f)^2 + a*e*f^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + 2*a*e*f^2*tan(d*x)^2*tan(1/2*c)^2*tan(1/2*d
*e/f)^2 + 2*a*e*f^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + a*e*f^2*tan(1/2*c)^4*tan(1/2*d*e/f)^2 + a*e
*f^2*tan(d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*e*f^2*tan(1/2*d*x)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + 2*a*e*f
^2*tan(1/2*c)^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*f^3*x*tan(d*x)^2*tan(1/2*d*x)^2 + 2*a*f^3*x*tan(d*x)^2*tan(1
/2*c)^2 + 2*a*f^3*x*tan(1/2*d*x)^2*tan(1/2*c)^2 + a*f^3*x*tan(1/2*c)^4 + a*f^3*x*tan(d*x)^2*tan(d*e/f)^2 + a*f
^3*x*tan(1/2*d*x)^2*tan(d*e/f)^2 + 2*a*f^3*x*tan(1/2*c)^2*tan(d*e/f)^2 + a*f^3*x*tan(d*x)^2*tan(1/2*d*e/f)^2 +
 a*f^3*x*tan(1/2*d*x)^2*tan(1/2*d*e/f)^2 + 2*a*f^3*x*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + a*f^3*x*tan(d*e/f)^2*tan(
1/2*d*e/f)^2 + a*e*f^2*tan(d*x)^2*tan(1/2*d*x)^2 + 2*a*e*f^2*tan(d*x)^2*tan(1/2*c)^2 + 2*a*e*f^2*tan(1/2*d*x)^
2*tan(1/2*c)^2 + a*e*f^2*tan(1/2*c)^4 + a*e*f^2*tan(d*x)^2*tan(d*e/f)^2 + a*e*f^2*tan(1/2*d*x)^2*tan(d*e/f)^2
+ 2*a*e*f^2*tan(1/2*c)^2*tan(d*e/f)^2 + a*e*f^2*tan(d*x)^2*tan(1/2*d*e/f)^2 + a*e*f^2*tan(1/2*d*x)^2*tan(1/2*d
*e/f)^2 + 2*a*e*f^2*tan(1/2*c)^2*tan(1/2*d*e/f)^2 + a*e*f^2*tan(d*e/f)^2*tan(1/2*d*e/f)^2 + a*f^3*x*tan(d*x)^2
 + a*f^3*x*tan(1/2*d*x)^2 + 2*a*f^3*x*tan(1/2*c)^2 + a*f^3*x*tan(d*e/f)^2 + a*f^3*x*tan(1/2*d*e/f)^2 + a*e*f^2
*tan(d*x)^2 + a*e*f^2*tan(1/2*d*x)^2 + 2*a*e*f^2*tan(1/2*c)^2 + a*e*f^2*tan(d*e/f)^2 + a*e*f^2*tan(1/2*d*e/f)^
2 + a*f^3*x + a*e*f^2)

Mupad [F(-1)]

Timed out. \[ \int \frac {\cos ^3(c+d x)}{(e+f x)^2 (a+a \sin (c+d x))} \, dx=\int \frac {{\cos \left (c+d\,x\right )}^3}{{\left (e+f\,x\right )}^2\,\left (a+a\,\sin \left (c+d\,x\right )\right )} \,d x \]

[In]

int(cos(c + d*x)^3/((e + f*x)^2*(a + a*sin(c + d*x))),x)

[Out]

int(cos(c + d*x)^3/((e + f*x)^2*(a + a*sin(c + d*x))), x)